Two-time scale subordination in physical 
processes with long-term memory 



Aleksander Stanislavsky , Karina Weron ^ 

^Institute of Radio Astronomy, Ukrainian National Academy of Sciences, 
4 Chervonopraporna St., 61002 Kharkov, Ukraine 

^Institute of Physics, Wroclaw University of Technology, 
Wyh. Wyspiankiego 27, 50-370 Wroclaw, Poland 



Abstract 

We use the two-time scale subordination in order to describe dynamical processes in 
continuous media with a long-term memory. Our consideration touches two physical 
examples in detail. First we study a temporal evolution of the species concentration 
for the trapping reaction in which a diffusing reactant is surrounded by a sea of 
randomly moving traps. The analysis is based on the random- variable formalism 
of anomalous diffusive processes. We find that the empirical trapping-reaction law, 
according to which the reactant concentration decreases in time as a product of an 
exponential and a stretched exponential function, can be explained by the two-time 
scale subordination of random processes. Another example is connected with a state 
equation for continuous media with memory. If the pressure and the density of a 
medium are subordinated in two different random processes, then the ordinary state 
equation becomes fractional with two time scales. This allows one to arrive at the 
state equation of Bagley-Torvik type. 

Key words: Subordination, Nonexponential relaxation. Trapping reaction. 
Anomalous diffusion, Continuous medium with memory 
PACS: 02.50.-r, 82.40.-g, 02.50.Ey 



1 Introduction 



The problem of trapping reactions lias a long and active history tracked in 
the literature (see, for example, p!|2|3] ). The reaction dynamics has been also 
studied over the past decade [H[5p|7|5] . In the traditional version a reactant 
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(A) walks in a medium randomly doped with static traps (B) so that when 
they meet, the reactant disappears {A + B ^ B). This has served as a start- 
ing point for the formulation of several models describing the behavior of 
more complex systems. An important variation of the basic trapping prob- 
lem is that the traps become diffusive too. When both species [A and B) are 
(sub)diffusive, the temporal evolution of the system looks like a random mo- 
tion (of walker A) in a random sea of traps (trapping random walkers B). One 
approach for the analysis of such processes is based on the continuous time 
random walk (CTRW) in which space random jumps follow among subsequent 
random waiting times PITU] . Another approach is stated on the fractional dif- 
fusion equation, which describes the probability density of finding the particle 
at position x at time t [TTp2] . Recently, the relationship between the CTRW, 
stable distributions and the fractional calculus has been established exactly 
[TBUl]. It turns out that if the CTRW is represented as a subordination of 
a space random (Brownian or Levy) process by a inverse-time a-stable Levy 
process, the probability density of a walker position is governed by a fractional 
Fokker-Planck equation. The approach was also extended on anomalous re- 
laxation processes [T5|16|17|IT8] . The aim of the given paper is to present a 
development of the subordination formalism to the complex relaxation pro- 
cesses. We are going to consider a time evolution of the concentration in the 
trapping reactions when the relaxation function takes the form of a product of 
an exponential and a stretched exponential function. The experiment support- 
ing this empirical law is described in the paper of Djordjevic [20]. The author 
studied the problem of an electron transfer from the methyl viologen radical 
cations to the colloidal platinum particle in a water solution. In that reac- 
tion the monocations MV"*" became the dications MV^"*" and the dependence 
of MV"*" concentration on time appeared to follow the combined exponential- 
stretched-exponential decay law. In particular the experiment has stated the 
relaxation curve following logC(t) = -0.0011825t - 0.0608^° *^ such that the 
cautious estimation of the exponent of the second term is 0.56 ± 0.06. Next, 
Djordjevic also noticed that the contribution of the second term is about 90% 
of log C{t). One of the aims of his paper was to verify experimentally the 
fact of factorization. The task was realized with success. We believe that this 
experiment is good, and the data fitting gives the physical dependence clearly. 

The paper is organized as follows. In Sec. [2] we present important features of 
subordinate random processes with different subordinators. They are directly 
connected with the anomalous diffusion. In the dependence of the subordinator 
form the corresponding subordinate processes have a different evolution in 
time. In our analysis we use the peculiarity of relaxation functions. From the 
paper [IS] it is known what subordinator leads to the stretched exponential 
law. In Sec. [3] we show that for the relaxation law to take the combined form, 
the subordinator should contain two internal-time scales. One of them gives 
the exponential relaxation, and another contributes the stretched exponential 
function. In the spirit of the consideration. Section H] is devoted to the study of 
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continuous media with a long-term memory. We consider such media that are 
characterized by the two time scales in the state equation. Section [5] presents 
a summary of results. 



2 Stretched exponential response and its subordinator 

If a complex physical system consisting of identical objects (dipoles, charges 
and so on) undergoes an irreversible transition, say, from state A to state B 
at random instances of time, then the transition can be characterized by the 
probability in the form 



where the non-negative quantity r{y) denotes the system's transition rate, i.e., 
the transition probability intensity for transition of the system as a whole. This 
formula simply follows from a two-stage master equation (see, for example, 
|21]). The probability Pr(^ > t) shows that a considered object will remain 
in state A until time t . In general, the transition rate r{y) is time-dependent 
because of random impacts affecting each object. If one knows the explicit form 
of r{y), the value Pt{6 > t) can be derived. In the simplest case, when the 
quantity r{y) = Bq = const is time-independent, the above formula recovers 
the classical exponential relaxation law 



where oOp = bo is a characteristic material constant. The exponential law, 
however hardly ever found in nature, is widely accepted for description of 
various relaxation data. Such a model of relaxation assumes that the relaxation 
rate (inverse of the relaxation time) of each object is constant. Although it 
may be different from one object to others, the mean of the effective relaxation 
rate (representing the inherent stochastic nature of the relaxation process) has 
to take a finite value. Nevertheless, the extensive experimental investigations 
in a wide frequency-time domain have shown P^E3] relatively large deviations 
from the exponential relaxation law. It has been found that for many materials 
the relaxation response follows the stretched exponential pattern 

0(t) = e-(^*'*)", 0<a<l. 

In this case, on contrary, the system's transition rate is essentially time de- 
pendent, and the mean of the effective relaxation rate becomes infinite 




Pr(0 >t) = e 
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The time dependence of the system's transition rate as well as the statistical 
properties of the effective relaxation rate clearly depend on the defect-diffusion 
relaxation mechanism in the system under consideration pjj. Hence, the ex- 
ponential and the nonexponential relaxation can be modeled by means of the 
diffusive behavior of the systems as a whole. Of course, the exponential relax- 
ation and nonexponential one correspond to different diffusion processes. 

Consider a sequence Tj, i = 1, 2, . . . of non-negative, independent, identically 
distributed random variables which represent waiting-time intervals between 
subsequent jumps of a particle. The random time interval of n jumps in space 
is written as 

n 

Tin)=J2T,, T(0) = 0. (1) 

i=l 



The number of the particle jumps performed up to time t > is determined 
by the largest index n for which the sum of n interjump time intervals does 
not exceed the observation time t 



Nt = max{n : T{n) < t}. 



The process A^^ is called the renewal process, or else, the counting process. 
The position of the particle (i.e., the total distance) after Nt jumps becomes 

Nt 

R{Nt)=J2R,, (2) 

i=l 



where Ri are independent, identically distributed variables giving both the 
length and the direction of the z-th jump. The process ([2]) is just known as 
the CTRW. 

Assume that the time intervals Tj belong to the domain of attraction of a 
completely asymmetric stable distribution with the index < a < 1. The 
generalization of the central limit theorem yields the continuous limit of the 
random sum ([T]), namely 

a-'^/''T{[aT]) 4 U{t) as a ^ oo, 

where U{t) is a strictly increasing a-stable Levy process, a > parameter, [x] 

denotes the integer part of x and "-4" means "tends in distribution" . Similarly, 
let the jumps Ri belong to the domain of attraction of a 7-stable distribution 

Sy,p{x), < 7 < 2, < 1 SO that 

a~^^'^ R{[aT]) A X{t) as a 00, 
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where X{t) is a 7-stable Levy process. If 7 = 2, the latter process is the 
classical Brownian motion. Both the process U{t) and the process X{t) are 
indexed by the internal time r. The time is not the real, physical time. In order 
to find a particle position at the observable time t, we have to introduce the 
notion of the inverse-time a-stable Levy subordinator Vt relating the internal 
T and the observable t times 

a-^'Nat AVt = inf{r : f/(r) > t} as a ^ 00. (3) 

Then the continuous limit of the CTRW process obtains the following form 

a-^h^^Nat) ^ (a")-i/^i?([a°\/t]) 4 X{Vt) as a ^ 00 (4) 



known as the anomalous diffusion process. In other words, the scaling limit of 
the CTWR leads to the anomalous diffusion process X{Vt) in which the parent 
process X{t), replacing the random discrete-time jumps R{n), is subordinated 
by the directing process Vt, replacing the counting process Nt. It has been 
rigorously proved that the probability density of X{Vt) is the solution of well- 
known fractional diffusion equation [TSUI]. As the processes X{t) and Vt are 
independent, following the total probability formula, the probability density 
Pa{x,t) of X{Vt) can be written as 

00 

Pa{x,t) = / fa{T,t) g{x,T)dT , 







where fa{r,t) and g{x,T) are the probability density of Vt and X{t) respec- 
tively. Similarly, the Fourier transform pa{k,t) = {exp{ikX {Vt))) and the 
Laplace transform pa{k,t) = (exp(— fcX(Vf))) are given by 



Pa{k, t) = J fair, t) g{k, t) dr , (5) 


00 

Pa{k,t) = / fa{r,t)g{k,T)dT. (6) 



Here k > has the physical meaning of a wave number. 

The analysis of the properties of the diffusion front (i.e., the asymptotic dis- 
tribution of the particle position in time t) allows one to determine [TBPT] 
the one-parameter Mittag-Leffler relaxation function (corresponding to the 
Cole-Cole law in the frequency domain), namely 

cp{t) = (e-'^^^^*)) = E,(-c«,,Ft"), (7) 



5 



where Ea{x) = J2kLo^"' /^{^(^ + 1) is the Mittag-Leffler function, Cq,^ is con- 
stant (cUp = {ca^'jk'^Y^"). For the relaxation function to take another form, the 
subordinator of the process X(t) should be changed. 

With that end in view the papers P^fTS] show that if the subordinator is a 
fully asymmetric Levy a-stable process with the self-similar property 

^/°) £ (8) 

where "=" reads "equal in law" , the type of relaxation remains simple expo- 
nential 

(^e-'''^''''^ = e-^"'"\ 0<a<l. (9) 



This means that there are different scenarios leading to the exponential term 
in the relaxation function. It can be the Brownian as well as a deterministic 
process. 

Nevertheless, if a new subordinator v["^ is defined as [18] 

Fj^^^cty/"), oO, (10) 



then the relaxation law decays otherwise. This process is constructed on the a 
stable Levy process Vt""^ multiplying both sides of its self-similarity property 

V},"^ I {ctY/»vl-\ OO, (11) 



by (ct)^"^/". Both processes and V^j" have different properties. The old 
one Vf°'^ is an 1/a self-similar Levy process, and this means that it is an a 
stable Levy process. While the new one is an a stable process, not Levy (i. 
e. it does not have independent increments). Note that any Levy process is 
a process with independent and stationary increments, whereas the a stable 
is a process for which all finite dimensional distributions are stable. Hence, 
the a stable Levy process v}"''^ has independent and stationary increments as 
well as all finite dimensional distributions stable. It is a broad class of random 
processes (see Table in [I8j). 

It would be useful to notice that the self- similarity of the process v}"'^ can be 
also written as 
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where the constant a is dimensionless. Then the corresponding property for 

— (a) 

the new process VI takes the form 
Kt = , a > 0. 



This once again confirms their different behavior. 

Now consider the subordinated process X{V^^^) in which for simphcity a ran- 
dom process X(t) is ordinary Brownian, and the directing process is v[°'\ Let 
the subordinated process X(v['^^) with such a subordinator have a probabihty 
density pa{x,t). The relationship between the probabihty densities of parent 
and directing processes is expressed in the integral form 

OD 

Pa{x,t) = j pi{x,t^)dga{0 : (12) 


where the probability distribution QaH) is described by the Laplace transform 

oo 

G^{s) = J exp{-s^}dg^{0 = exp{-(As)"} 





with s > and A > 0. The probability density pi{x,t) describes the normal 
diffusion. The anomalous diffusion gives the stretched exponential 

function of relaxation 

L-wr^X^^-c^f-k-^ 0<a<l, OO. (13) 



In the next section we will demonstrate how the combined laws of relaxation 
can appear in the evolution of relaxing physical systems due to a further 
development of the subordination approach. 



3 Subordination by two random processes 



Let two independent sequences {t/^^} and {T^^^^}, i — 1,2, . . . consist of non- 
negative, independent and identically distribTitcd random variables. The ran- 
dom variables Tj;^^ and t/^-* are attracted to a-stable laws with different in- 
dices. Let the temporal variables T^^\n) and T^'^\n) be a sum of the sequences 
of time intervals, t/^-* and t/^-* respectively. The counting process describes 
the number of the particle jumps performed up to time t. The position of the 
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particle after n jumps is a sum of the jumps The total distance reached 
by a particle during time t is defined by the number of the jumps by means 
of the counting process N^. In this case the process A''^ is something like the 
operational time. Since the sequences of the time intervals are independent 
on each other, their subordinators will be such too. To pass to the continuous 
limit, the probability density of the obtained subordinate process should be 
expressed as 



oo oo 

n + T2 



fa{r2, t) g(^x, ^ ^ ) ffsin, t) dTi dT2 , (14) 







where the variables ri and T2 correspond to the two-time scale subordination 
for the parent process with the probabihty density g{x, r). Here, we also have 
an anomalous diffusion. 

In the case of the nonbiased random walk the process X{t) belonging to the 
class of 7-stable Levy process has the following characteristic function 



If the parent process is directed by the two subordinators, then the relaxation 
function takes the form 



00 00 

-c^fcT(ri+r2)/2 







As a result, we obtain a combined process. 

Now, we look at the trapping reaction with the randomly moving traps. Ac- 
cording to the experimental data [20j, one can prepare the corresponding ex- 
perimental system so that after activating the traps, the reactant concentra- 
tion will relax in correspondence with a combined law. This law may be a 
product of some (in the simplest case, two) well-known relaxation laws (such 
as exponential, stretched exponential and so on). From the above considera- 
tion it follows that the relaxation function for the trapping reaction can be 
expressed in terms of 



= e-''"'''''-e-''^'^''\0<a,(3<l, 0^,0/3 >0. (16) 
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To sum up, the relaxation function is written as a simple product of the ordi- 
nary exponential and stretched exponential functions. Two random processes 
are present at the reaction: 1) reactant walks in a medium; 2) the traps appear 
randomly. If these processes are independent, then their contribution in the 
relaxation function can be a product of two relaxation dependencies. One of 
them gives a simple exponent, and another tends to a stretched exponential 
function. The indices a and /3 permit one to distinguish the random processes, 
but we do not define concretely what component (reactant or traps) leads to, 
for example, the stretched-exponent contribution because this depends on the 
experimental situation. We describe the most general picture of this phenom- 
ena from the probabilistic formalism of limit theorems. 

Why is subordination important here? Because the reactant walks is (sub) diffu- 
sive, and (sub) diffusion is a result of subordination of one random process by 
another [25] . In like manner this relates to traps. The nontriviality of our anal- 
ysis is due to the fact that the effective stochastic process in time at a fixed 
point in space becomes non-Markovian. For a non-Markovian process, calcu- 
lation of any history dependent quantity such as functional of trajectories 
is extremely hard barring a few special cases (see more details, for example, 
in [25]). Over past decade the probabilistic formalism of limit theorems has 
made a very great advance on the analysis of non-Markovian processes. We 
just apply it for the treatment of Djordjevic's experiment. 

It should be emphasized that the subordination approach answers an impor- 
tant question: what a stochastic mechanism stands behind the combined law of 
relaxation in the studied system. We believe that this is caused by a multi-time 
scale subordination of random physical processes in such systems. Moreover, 
this is not the only physical development of multi-time scale subordination. In 
the next we consider how the subordination can influence on a state equation 
of continuous media. 



4 State equation for continuous media with memory 

Any continuous medium at any point x = {x\^X2^xz) G and any time t 
is determined by the velocity vector v(x, t) = (fi,f2,f3) , the pressure p(x, t) 
and the density p(x, t). The characteristic quantities are connected with each 
others by means of the Navier-Stokes equation 

dv 

p— + Vp = /iAv, 
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the continuum equation 

— + div pv = 
at 



and the state relation 
Fip,p,v,yi,t) = 0, 



where V = ^) is the Hamiltonian operator, n the viscosity, A the 

Laplace operator, and F is an operator (see [27]). 

Following [28] , for a start we assume that the relaxation process in a continuous 
medium obeys a simple exponential law 

dh/dt = -{h - ho)/e or h = ho + {h{0) - ho)e-'^^, (17) 



where we do not define concretely the nature of relaxation, and 6 is the re- 
laxation time of the parameter h characterizing this medium. Next, a wave 
(sound or something like) is propagated in the medium. If the wave period is 
T ^ 9^ the media has time to accommodate oneself to wave changes, and the 
wave will advance with the velocity Cq = {dp/dp)ho (here ho corresponds to 
a value of h in this wave). If on the contrary T 6, then the parameter h 
will be "frozen" to its value /iqo without any wave, and the wave velocity is 
= (dp/dp) hoo- Suppose that the wave results in small perturbations of the 
density p = Po + p' and the pressure p = pq + p' for the given medium. 

Let the state relation take the form p = p{p, h). Expand it in a series about 
small perturbations: 



P = P{p, h)^P0+ {^)JP - Po) + {gj^)jh - ho). (1 

As the equilibrium state has ho = go{p), therefore 



pip, ho) ^Po+ (I), (P - Po) + {^)iho - hoo). (19) 



It follows from this that 

^dp^ ho ^dp^ hoo^ ^Qji) po Qp ^ ^ 
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Next, we substitute the expression (l20i) to the expansion ([T8|) and differentiate 
the obtained result with respect to t. Using the equation (fT7|) one can find 

dp' _ 2<^_(^\ h-hp 

dt ~ dt ydg)po 9 ■ ^ ^ 



The expression (12T|) shows a connection between increments of p' and p' in 
the wave and the deviation of the parameter h — hp from its equihbrium value. 
In order to get an equation depending only on and p', we divide Eq. (fT8|) on 
T and sum the result with Eq. fl^T]) . Thus we arrive at the sought-for equation: 

dp^ P^_ 2dp^ 2P^ 

dt^ e~^^ dt ^^'e ■ ^ ' 



It is not difficult to establish that this equation is equivalent to the integral 
relation 

P' = cy + {cl- cl) } exp (-^) ^ dt, . (23) 

— oo 



From Eq.f l23p it is seen that the medium has memory effects. 

For the derivation of Eq. (!23|) we accepted the relaxation equation of type 
( [T7|) . As a result, the kernel of Eq.( l23|) has the exponential form. Generally, 
the relaxation equation can be more complex (for example, for polymers it 
is non-exponential), and therefore the corresponding kernel may be other. 
In particular, for fractional systems it is written as a power function, and the 
corresponding state equation takes the fractional form of differential equations. 
Consider the feature in more details below. 

A broad class of continuous media with a long-term memory (for example, 
polymer ffuids, viscoelastic materials, etc. |29130J) is described more ade- 
quately by the state equation in the fractional form of type 

p^_ ,p^ 
dt- + ^1 " °° dt^ + ^ ' 



where d^/dt"' is the fractional derivative in the sense of Caputo [12]. Here Cq is 
the wave velocity, when the wave period is more than relaxation time 6i, and 
Coo is the wave velocity in the opposite case. To obtain Eq. (!2^ from Eq. (!22|) . 
we present the pressure and the density as a result of subordination, namely 
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Pa(x,t)=y" faiTi,t)p{x,Ti)dTi, (25) 


oo 

pp{-K,t) = j /^(r2,t)p(x,r2)rfr2, (26) 



where faijiji) (or fi3{T2,t)) is the probabihty density of the inverse-time a 
(or /3)-stable Levy subordinator with < < 1. In other words, both 
the pressure and the density of such a medium are subordinated in different 
ways. From the physical point of view the subordination accounts that there 
are temporal (random) intervals, when the pressure and the density does not 
change. In any other case the random jumps in density and pressure occur. In 
particular, if = 1, then the values change continuously. 



According to [12], the Laplace transform in time, for the pressure and the 
density in the medium with the long-term memory, gives 

p„(x, s) = s"-ip(x, s°), pp{^, s) = /-^(x, /). (27) 
Hence, Eq. (l22|) in the Laplace space reads 

p(x, s) 



sp(x, s) -pix, 0) + 



'1 



= clspi^, s) - clpi^, 0) + cl^^. (28) 
The subordination in different ways means that we should accept 



7 ( gp(x,ri) p{^,Ti) \ 





oo 



2 <9p(x,r2) 2P(x,1-2) 







cL ' + ^P-Jf^ ) M^2,t) dr,. (29) 



Then the Laplace transform of the latter expression takes the form 

."-i(."p(x,s")-Kx,0)) + ^^ = 
cls^-\s^p{^. /) - P(x, 0)) + clP^. (30) 
After simple algebraic transformations we have 
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s°p„(x,s)-s"-V(x,0) + 



Pa IX, s) 



(31) 



and the Laplace inverse of this expression leads directly to Eq.( l24|) . The sub- 
ordination of type fl25|) and ( 126|1 does not change the form of the Navier-Stokes 
equation, but in the continuum equation the time derivative for the density 
becomes fractional. 

Finally, we consider the situation where both the Navier-Stokes equation and 
the continuum equation are one-dimensional, namely 



dv dpa d^v 

P'm^^ = ^d^^ ^^'^ 



d^ps dv 

'W ^ " °' ^ ^ ^ °' ^^^^ 

For a = (3 and cq << 1 Eq.( l24l) tends to the following expression 

^i^+P^ = ^i^^^- (34) 



Let the viscosity /i be equal to zero. From Eqs. (l32|) and ( 133|) we have 
dt at" ~ 'dx^' 



Using the state equation flM|) . the equation for the pressure takes the form 

^^^'m^-^''-^- ^^^^ 



The presented case is interesting because it leads to the fractional wave- 
diffusion equation ( l35ll . although the Navier-Stokes equation is not fractional. 
The fractional derivative appears in Eq. fl35|) . due to the fractional state equa- 
tion and the fractional continuum equation. 



5 Conclusions 



We have shown how the combined law of relaxation can be derived from the 
diffusion model based on the CTRW consideration and illustrated the role of 
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the two-time scale subordination. Our analysis demonstrates that the type 
of relaxation function depends entirely on the subordinator. The subordi- 
bators may be different and enough complicated. If the operational time is 
an inverse-time Levy ct-stable subordinator, the relaxation function becomes 
Mittag-Leffler's. The operational time of type VI "transforms" the relaxation 
function into the stretched exponential form. The multi-scale subordinator al- 
lows one to combine some laws of relaxation as their product. The result 
supports the conclusion that the subordination as a transformation from the 
physical time t to the operational time r is responsible for the anomalous 
properties of complex systems. Using this approach to a similar problem in 
continuous mechanics, we notice a special role of the state relation. Due to the 
multi-scale subordination the relation takes the fractional form. Consequently, 
the continuum equation becomes mathematically more complicated than one 
in the ordinary mechanics. 
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